Interband tunneling in nanowires with diamond cubic or zincblende crystalline 

structure based on atomistic modeling 
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We present an investigation of band-to-band tunneling in nanowires with a diamond cubic or 
zincblende crystalline structure, as a function of cross section and longitudinal electric field. Results 
are obtained from quantum transport simulations based on non-equilibrium Green's functions with a 
tight-binding atomistic Hamiltonian. Interband tunneling is extremely sensitive to the longitudinal 
electric field, to the nanowire cross section, through the gap, and to the material. We have derived 
an approximate analytical expression for the transmission probability based on WKB theory and on 
a proper choice of the effective interband tunneling mass, which shows good agreement with results 
from atomistic quantum simulation. 
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Conventional field-effect transistors (FETs) are based 
on the modulation of thermionic injection of charge car- 
riers from the source by tuning the potential barrier be- 
tween channel and source via the voltage applied on the 
gate electrode. A crucial parameter in FETs is the sub- 
treshold swing S - defined as the gate voltage variation 
required for a tenfold increase of the drain current. In 
the ideal case of perfect electrostatic control, when the 
potential in the channel exactly follows the gate voltage, 
S is limited by the tail of the Maxwell-Boltzmann distri- 
bution to fcsTlnlO, where fee is Boltzmann's constant 
and T is the temperature. 

At room temperature, the intrinsic limitation on S is 
therefore S > 60 mV/decade, which represents a real 
threat to continuous integration of semiconductor tech- 
nology, due to conflicting constrains. Indeed, supply volt- 
age (Vdd) reduction must accompany device size scaling 
to keep power consumption per unit area (oc V^ D ) un- 
der control, and digital logic still requires the ratio of the 
current of the device in the on state to be at least 10 4 
larger than than in the off state, effectively requiring the 
threshold voltage V t h > 45. 

A promising option to reformulate this trade-off is to 
use a device based on a modified operating principle such 
as the tunnel FET (TFET) Q. In a TFET the gate 
voltage allows to modulate the band-to-band tunneling 
(BTBT) current between source and channel, by modu- 
lating the width of the energy window in which tunneling 
can take place. In this way, both ends of the allowed en- 
ergy window for transport are defined by conduction and 
valence band edges, and the previously discussed limita- 
tion on S is therefore removed, allowing further reduc- 
tion of Vd d ■ Several experimental and theoretical works 
have been dedicated to explore various options for tunnel 
FETs. @-|| 

Two aspects are important for good performance of 
TFETs. On the one hand, currents in the on state must 
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be large for obtaining small delay times, shifting one's 
preference to lower gap semiconductors, when BTBT can 
be large. On the other hand, electrostatic control of the 
potential in the channel via the gate voltage must be 
ideal, not to unduly lose the advantage on S, as can be 
obtained in ultra-thin-body or gate-all-around devices. 

In this letter, we investigate the dependence of band- 
to-band tunneling on nanowire diameter, material, and 
longitudinal electric field in the interband tunneling bar- 
rier, in order to understand their role in determining the 
achievable BTBT current, to allow comparison of dif- 
ferent material systems and channel geometry, and to 
provide reliable quantities for use in device models at a 
higher level of abstraction. 

We perform atomistic simulations based on a sp 3 d 5 s* 
Hamiltonian for diamond cubic or zincblende crystals. 
In particular, we focus on germanium and indium ar- 
senide - which arc interesting candidates for TFETs for 
their relatively small gap as bulk materials - and on sili- 
con, as a reference material for maximum ease of integra- 
tion with CMOS technology. Nanowire transport is com- 
puted in the framework of non-equilibrium Green's func- 
tions (NEGFs). Finally, wc propose a simple analytical 
expression derived from the Wentzel-Kramcrs-Brillouin 
(WKB) approximation with an ad-hoc tunneling effec- 
tive mass, that is able to quantitatively reproduce re- 
sults from atomistic quantum transport simulations and 
can be used for device modeling at a higher level of ab- 
straction. 

In Fig. [TJa) we show an example of a transversal cross 
section of a zincblende nanowire with the axis oriented 
in the [001] direction. We use a tight-binding Hamilto- 
nian with a sp 3 d 5 s* first-nearest neighbour representa- 
tion, where 10 atomic orbitals are considered (spin orbit 
coupling is neglected). This representation is among the 
most accurate descriptions available and it is the most 
used for refined numerical calculations 041] ■ It also gives 
quantitative agreement with experimental results (lOJ. 

We have developed a tool to handle generic Hamiltoni- 
ans with an arbitrary number of atomic orbitals. In order 
to speed-up the calculation of the transmission proba- 
bilities, we have optimized the tool exploiting the four- 
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Figure 1: Example of the investigated structures and its main 
properties: InAs squared nanowire with diameter of 3.4 nm: 
(a) lattice structure of the (001) plane transversal to the 
tunneling direction of the wire; (b) band structure of the 
nanowire; (c) band edge profile considered for the calcula- 
tion of the tunneling coefficient; the tunneling energy window 
is highlighted; (d) transmission coefficient in the tunneling 
window as a function of energy for different values of the lon- 
gitudinal electric field. 



atomic-layer periodicity of zincblende materials in the 
[001] direction. 

The computation of transmission is based on the use of 
non-equilibrium Green's functions, using a recently de- 
veloped closed-form method for the calculation of the 
lead self-energies. [l2j Due to the complexity of the in- 
vestigated structures, a closed-form self-energy scheme 
was required to reduce the computing time with respect 
to iterative procedures. In particular, the adopted proce- 
dure has been demonstrated to be ten times faster than 
the Sancho-Rubio iterative method. [l3[ 

The transmission probability is obtained with the well 
known formula 



T = Tr[G r -r L -G a -r R ], 



(1) 



where G r ' a are the retarded and advanced Green's func- 
tion matrices and Tl^r represent the tunncling-rate ma- 
trices for the left and right lead, respectively. A detailed 
description of the NEGF transport theory can be found 
in the literature. fl4| 

To investigate the possible use of nanowires with 
zincblende crystalline structure as channels for tunnel 
FETs we have computed the transmission coefficient 
due to band-to-band tunneling. We first compute the 
dispersion relations of the nanowire under investigation 
as shown in Fig. H^b), in order to extract the energy 
gap. Then, we apply a linear potential as shown in 
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Figure 2: The plateau of transmission probability as a func- 
tion of the longitudinal electric field is shown for nanowires 
along the [001] direction, with square cross section of side 
I = a ^-N with N = 4 (a), N = 6 (b) , N = 8 (c), where a 
is the lattice constant (an = 5.43, 5.65, 6.05A) for Si, Ge, and 
InAs, respectively). In every plot the transmission coefficient 
is shown for three different materials: Si, Ge and InAs. (d) 
dependence of the energy gap on I, for the different materials, 
extracted from the dispersion relation. 



Fig. HJc) with and end-to-end potential drop larger than 
Eg in order to open an inter-band tunneling window of 
0.3 eV. This ensures limited mismatch and a transmis- 
sion probability almost independent of energy, as shown 
in Fig. [ljd), where trasmission probability exhibits a 
plateau in the interband tunneling energy window. In 
the following analysis, we use the value of the transmis- 
sion coefficient in such plateau as a metrics for comparing 
interband tunneling among different materials and as a 
function of the applied longitudinal field. 

Nanowire tunnel FETs pose different tradeoffs: a 
smaller nanowire cross section implies a better electro- 
static control, and therefore a better subthreshold behav- 
ior. On the other hand, it also implies a larger gap, and 
therefore a smaller tunneling current in the on state. To 
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explore this tradeoff from a quantitative point of view for 
different materials, we consider nanowires in which prop- 
agation occurs along the [001] direction, and the lateral 
surfaces are on {110} planes, as shown in Fig. Ufa). 

The side of the square cross section I is determined 
by an integer parameter N according to the formula 
I = cio^-N where ao is the lattice constant for the ma- 
terial under investigation (a = 5.43, 5.65, 6.05A) for Si, 
Ge and InAs, respectively). In Fig.[5]we show results for 
N = 4,6, 8, corresponding to a side I between roughly 
1.5 nm and 3 nm depending on the material. As can 
be seen in Fig. [2jd), with increasing cross section area 
the gap rapidly decreases, causing a very steep increase 
of the tunneling probability. In addition, as can be seen 
from Figs. [2Ia)-(c) the transmission coefficient increases 
by several orders of magnitude for an increase of the lon- 
gitudinal electric field of just one order of magnitude. 
In the case of InAs, in particular, for the nanowire with 
I w 3 nm, the transmission probability can be very high, 
opening the possibility for very large ON currents. 

As we pointed out before, full band calculations, as the 
ones presented here, are computationally very demanding 
for diamond cubic and zincblende nanowires. It is there- 
fore interesting to verify whether simplified expressions 
are reliable in predicting the behavior of the tunneling 
coefficient. In Fig. [3] a comparison between the trans- 
mission probability plateau obtained from the atomistic 
NEGF calculation and the WKB results is shown. The 
WKB transmission coefficient is obtained with the for- 
mula 



Iwkb = exp 



2m 



Fzdz 



(2) 



where F is the longitudinal electric field and m = mo ■ m» 
is the effective mass in the interband tunneling barrier 
plotted in Fig. [3Jd) (mo is the electron mass at rest). 
The potential profile is linear with the slope given by 
the applied electric field, as shown in Fig. [ljc). The 
integration region is determined by zq = E g /F. After 
few manipulations the formula @ can be approximated 
as follows: 
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Figure 3: Comparison between the transmission probability 
plateau obtained from the atomistic NEGF calculation (solid 
lines) and the WKB results (dotted lines). Nanowires with 
three different cross section sides / are considered: I = ao^N 
with N = 4 (a), N = 6 (b), N = 8 (c). The values of effective 
masses used for the WKB approximation are shown in (d); 
they have been extracted from the energy-band profile of the 
corresponding atomistic calculation for InAs and as a fitting 
parameter for Ge and Si. 
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20 lE, 
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where E g is expressed in eV and F in V/nm. With (|3|) we 
fit the numerical results for InAs and Ge, choosing the 
effective masses in two different ways: for InAs, which 
exhibits two clear well-separated conduction and valence 
bands, the effective mass is simply obtained from the 
second derivative of the first conduction band profile. 
For Ge the effective mass is just a fitting parameter, be- 
cause multiple quasi-degenerate conduction and valence 
bands are present. The values of the masses are shown 
in Fig. Eld). For Si, in order to reproduce the atomistic 
NEGF results for the tunneling probability, we have mod- 
ified the functional dependence of the WKB transmission 



probability on the electric field using the formula 



F 



WKB 



exp 



20 



E}, 



p0.71 



We have then used the effective mass as a fitting pa- 
rameter (values shown in Fig. Eld)). As can be seen, 
the simplified formulas of Eqs. (|3I4[) reproduce quite ac- 
curately full-band calculations for InAs, Ge and Si re- 
spectively. A more complex implementation of the WKB 
approximation could in principle be used to obtained 
better accuracy but our results show that the sim- 
plest version can still provide a very good quantitative 
agreement. Therefore, first-principle numerically accu- 
rate calculations can be reproduced by simple analytical 
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formulas, with effective masses extracted from the band 
profile of the atomistic structure or from fitting proce- 
dures. This makes self-consistent calculations of currents 
a much less demanding computational task compared to 
full-band NEGF calculations. 

We have investigated diamond cubic or zincblende 
nanowircs for different materials and for different cross 
sections up to a side of 3.4 nm using atomistic Hamilto- 
nians based on sp 3 d 5 s* tight-binding representation. 

The extreme sensitivity of the tunneling probability to 
the electric field and to the cross section side provides 
a hint of the sensitivity of the tunnel FET ON current 
to the variability of the nanowire cross section and to 
the presence of random charged defects in the nanowire 
region where interband tunneling takes place. We have 
shown that tunneling probabilities between 0.01 and 0.1 - 



suitable for achieving reasonable ON currents in TFETs 
- can be obtained with longitudinal electric field close to 
1 V/nm (still smaller than the dielectric breakdown field) 
in InAs for I > 1.5 nm and for Ge for I > 3 nm. We have 
also shown that full-band calculations are quantitatively 
reproduced by simple formulas derived from the WKB 
approximation, considering an ad hoc interband tunnel- 
ing mass. 

Finally, the presented results - obtained with rigor- 
ous quantum transport simulations based on an atom- 
istic tight binding Hamiltonian - can be used as input 
for semiclassical modeling of electron devices, where in- 
terband tunneling has to be represented with a synthetic 
scattering rate. 
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